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Gamma-ray astronomy will play a crucial role in the investigation of nonthermal processes in 
the large scale structure of the universe. Particularly, galaxy clusters (GC) observations at this 
photon energy will help us understand the the origin of radio emitting high energy particles, 
the possible level of cosmic-ray (CR) pressure in intracluster environment, and the strength of 
intracluster magnetic fields. In addition here we point out the importance of these observations 
for a possible detection of cosmic shocks through 7-ray emission and for constraining their CR 
acceleration efficiency. We model spatial and spectral properties of nonthermal 7-ray emission 
due to shock accelerated CRs in GC and emphasize the importance of imaging capability of 
upcoming 7-ray facilities for a correct interpretation of any observational results. 

1 Introduction 

Cosmic shocks emerge during structure formation in the universe, as gravitationally driven gas 
infall onto collapsing objects becomes supersonic. The dissipation of kinetic energy at these 
shocks raises the temperature of the intergalactic gas. In galaxy clusters (GC), characterized by 
accretion velocities of order ~ 10 3 km s _1 , the gas temperature reaches ~ 10 8 K. Also, according 
to recent numerical simulations most of the low redshift baryons, observed at high redshift as 
Ly-a absorbers have been shock heated to a temperature of 10 5 — 10 7 K0. 

Astrophysical shocks are collisionless and it is believed that as part of the dissipation mech- 
anism supra-thermal populations of particles (cosmic-rays, CR hereafter) are generated via first 
order Fermi mechanism! Both CR protons and electrons can be accelerated. CR electrons 
are characterized by rather short radiative lifetime. For a typical intra-cluster environment, the 
main loss mechanism is inverse Compton (IC) for energies above ~ 150 MeV and Coulomb losses 
below that. Low energy, sub-relativistic CR protons also suffer Coulomb losses as electrons of 
the same energy. However, for relativistic protons the dominant energy loss mechanism up to 
the energy threshold for photo-pair production is p-p inelastic collisions with a timescale longer 
than a Hubble time. Therefore, once accelerated CR protons up to ~ 10 15 eV accumulate within 
large scale structures where they can be confined by fiG strong turbulent magnetic fields a. 

It is of interest to investigate the acceleration of CRs at cosmic shocks for a number of rea- 
sons. Firstly, GCs exhibit non-thermal radiation. This mainly consists of diffuse radio emission 
extending on Mpc scales. It is thought to be synchrotron radiation but there is no consensus as 
to the origin of the emitting relativistic electrons. In addition, excess of radiation with respect 
to thermal emission has been reported at both hard X-rayu and, although more controversially 
for the extreme ultraviolet part of the spectrum i'i. Secondly, if shocks acceleration operates 
efficiently, the proton component could bear a significant fraction of the total gas pressure, af- 
fecting the dynamics of cosmic structures I. Finally, CR electrons accelerated at inter-galactic 
shocks could contribute a significant fraction of the cosmic 7-ray background (CGB) @. 



In this contribution we address the role of 7-ray observations for investigating CR accel- 
eration at large scale structure shocks. The results presented here are based on a numerical 
simulation which is briefly described in In §|3] we address the contribution of cosmological 
CRs to the CGB and, based on EGRET experimental results, we constraint the efficiency of 
shock acceleration of CR electrons. In ^ we present the non-thermal spectrum at 7-ray energies 
for a Coma cluster prototype and in §|| we discuss the importance of 7-ray observations for 
investigating physical conditions in intra-cluster medium (ICM) environment. 



2 Numerical set-up 

We perform numerical simulations that model simultaneously the formation of the large scale 
structure and the acceleration, transport and energy losses of three CR species: namely shock 
accelerated (primary) protons and electrons and secondary e 1 * 1 generated in p-p inelastic collisions 
of the simulated CR protons and the thermal nuclei of the IGM. The simulation is fully described 
elsewhere! and here we shall only summarize the most important aspects of it. 

As for the large scale structure we adopt a canonical, flat ACDM model with total mass 
density O m = 0.3, vacuum energy density = 1 — O m = 0.7, normalized Hubble constant 
h = #o/100 km s _1 Mpc -1 = 0.67, baryonic mass density, f2& = 0.04, and spectral index and 
cluster normalization for initial density perturbations, n s = 1 and erg = 0.9 respectively. 

The CR dynamics, including injection at shocks, acceleration, energy losses/gains and spatial 
transport, is computed numerically through the code COSMOCRE^'B. CR protons are injected 
at shocks according to thermal leakage prescription. With the adopted parameters, a fraction 
about 10 -4 of the particles crossing the shock is injected as CRs. For shocks with Mach number 
between 4 and 10, which are responsible for most of the heating of the IGmH, this roughly 
corresponds to converting 30% of the shock ram pressure into CR proton pressure. The injection 
of primary CR electrons in the acceleration mechanism is here simply modeled by assuming that 
their number ratio to protons at relativistic energies be given by a parameter R e / p . Observations 
0111 typically indicate R e / p ~ 10~ 2 . In accord to the test particle limit of diffusive shock 
acceleration theory 13, the injected particles are redistributed in momentum as a power-law with 
slope related to the shock Mach number, M, as: q = 4/(1 — M~ 2 ). Each CR population is 
passively advected with the flow. Finally, energy losses/gains of the accelerated CRs are followed 
by solving numerically a "kinetic" equation written in comoving coordinates and integrated over 
a discrete set of logarithmically equidistant momentum intervals that define the numerical grid 
in momentum spaced. We follow CR protons in the momentum range between 0.1 GeV/c and 
10 6 GeV/c, and CR electrons and e ± between 15 MeV/c and 20 TeV/c. 



3 Contribution to the Cosmic Gamma-ray Background 

In Fig. gi (left) we illustrate the contribution to the CGB from CRs accelerated at cosmic 
shocks, together with observational data taken from the EGRET experiment &3( solid dots). The 
plot shows as a function of photon energy, e, the quantity 
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where j(s,z) is the computational-box-averaged spectral emissivity in units 'photons cm~ 3 s — 1 ' 
computed at each simulation red-shift, z, and at the appropriately blue-shifted photon energy 
e(l + z). In addition, r 77 is an attenuation factor due to photo-pair creation, 77 — ► e^, c 
indicates the speed of light and z max an upper limit of integration. We consider the following 
emission processes: IC emission of CR electrons scattering off cosmic microwave background 
photons (dot line), decay of neutral pions (it — > 77) produced in p-p inelastic collisions (dash 
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Figure 1: Left: Total 7-ray background flux produced by cosmological CRs (solid line), contributions from IC 
emission from primary CR electrons (dot line), 7r°-decay (dash line) and IC from secondary e (dash-dot line), 
and EGRET experimental data with their error-bars (solid dots) 13 Right: The IC 7-ray photon luminosity above 
100 MeV from individual clusters as a function of the cluster X-ray core temperature, T x . The flux is produced 
by CR electrons accelerated at accretion shocks onto GC. EGRET experimental upper limits are from Reimer 14 . 



line) and IC emission from secondary (dot-dash line). In fig. |l] the total flux (solid line) 
corresponds roughly to a constant value at the level of 0.2 keV cm -2 s _1 sr _1 throughout the 
spectrum. It is dominated by IC emission from primary electrons. A fraction of order 30% and 
10 % is produced by 7r°-decay and IC emission from secondary respectively. 

Remarkably, all three components produce the same flat spectrum, similar in shape to the 
observed one. This result is a reflection of the fact that the CRs distributions producing the 7-ray 
radiation were generated in strong shocks with a log-slope q ~ 40. However, the computed flux is 
only ~ 15 % of the observed CGB. It is difficult to imagine a higher contribution from 7r°-decay 
(dash line) and IC emission from (dot dash line). In fact, if more CR protons were produced 
at shocks, CR-induced shock modifications would actually end up reducing the population of 
7-ray emitting protons (and ). It is possible, however, at least in principle that CR electrons 
would be accelerated more efficiently than assumed here. In fact, by setting R e / p = 10~ 2 we 
assume that for strong shocks a fraction ij~4x 10~ 3 [log(p max /m e c)/log(4 x 10 7 )] (p m ax is the 
maximum momentum of the accelerated electrons) of the shock ram pressure is converted into 
CR electrons. This is an important matter because if one assumes that 77 is of order of a few 
percent then the origin of the unexplained component of the CGB could be accounted forH. 

However, the electron acceleration efficiency can be constrained by comparing the simulated 
clusters' p^-ray photon luminosity above 100 MeV to the upper limits set by the EGRET exper- 
iment Ej'O for nearby GCs. This is done in fig. [lp (right panel). The best-fit to the plotted 
data (solid curve) was found through a least-x 2 analysis to be: 

L 7 (> 100 MeV) = 8.7 x 10 43 (j^) f^)'' 6 ph s~\ (2) 

According to our calculations, the EGRET upper limits (particularly those relative to Coma and 
Virgo clusters) require that 77 < 0.8%. For the adopted proton injection efficiency this roughly 
translates into R e / P < 2 x 10 -2 . In any case, this implies an upper limit on the computed 7-ray 
flux of about 0.35 keV cm~ 2 s _1 sr _1 . Thus, according to our computation, cosmological CRs 
could contribute an important fraction of order ~ 25 % of the CGB. 



4 7 -ra y f° r a Coma-cluster prototype: synthetic emission maps and spectrum 



In this section we explore the spectral and spatial properties of non-thermal 7-ray radiation 
between 10 keV and 10 TeV due to shock accelerated CRs in GC. Clusters observations in 
this photon energy range are important for a number of issues which we address in detail in 
the discussion section. Among these are the following: correct measurements of intra-cluster 
magnetic field, estimate of CR proton pressure component, possible detection of cosmic shocks, 
estimate of CR acceleration efficiency. However, for a correct interpretation of the observational 
results it is important that the emission produced by different particle components be properly 
identified. In this respect, the results presented in the following are aimed at helping a correct 
diagnostic of future 7-ray observations. 

As non-thermal bremsstrahlung turns out unimportantly, we consider the following emission 
processes: 7r°-decay and IC emission from both primary CR electrons and secondary . We 
first compute the emission spectrum for the largest virialized object in the computational box 
which has an X-ray core temperature T x ~ 4 keV. The various emission components are then 
renormalized to the case of a Coma-like cluster. In particular, IC emission from electrons is 
rescaled according to equation ||. The total number of is rescaled assuming that these 
particles are responsible the for synchrotron emissionof Coma radio halo. For the purpose 
we took a radio flux at 1.4 GHz Si.iGHz = 640 mJyO and a volume average magnetic field 
(B) ~ 3/xG.ll! For the ensuing discussion it is not necessary to assume that the CRs responsible 
for the synchrotron emission are secondary . We know that these high energy particles exist 
from the observed radio flux and the morphology of radio halos suggest that they reside in 
the inner region of clusters. However, assuming that these are secondary e ± allows us to also 
normalize the rate of p-p inelastic collisions that determine the total 7r°-decay flux. 

The simulation results indicate that the radiation emitted by primary CR electrons and both 
7r°-decay as well as secondary originates in spatially separate regions. This is a reflection 
of both the spatial distribution of the emitting particles as well as the nature of the emission 
process. Thus, on the one hand, both and tt° are produced at the highest rate in the densest 
regions where both the parent CR ions and target ICM nuclei are most numerous. Consequently, 
IC from and 7r°-decay emissivities are strongest in the cluster inner regions and quickly fade 
toward its outskirts. On the other hand, because of severe energy losses, 7-ray emitting high 
energy primary electrons are only found in the vicinity of strong shocks where they are being 
accelerated. Thus, because the strongest shocks are located at the cluster outskirts rather than 
at its core where the ICM temperature is hightia, the spatial distribution of the emissivity is now 
reversed with respect to the previous case. In the four left panels of fig. || we present synthetic 
maps of the integrated photon flux above 100 keV (top panels) and 100 MeV (bottom panels). 
Among these four panels, the right panels are associated with IC emission from primary electrons 
and the left panels to IC from secondary (top-left) and 7r°-decay (bottom-left). As anticipated 
the emission from 7r°-decay (bottom left) and (top left) is confined to the cluster core. There 
it creates a diffuse halo which rapidly fades with distance from the center. On the other hand, 
IC emission from primary e~ is distributed over a much more extended area. Moreover, it 
is characterized by a strikingly rich and irregular structure. This is a direct reflection of the 
complex "web" of shocks that reside at the outskirts of galaxy clusters E3. 

After generating multifrequency synthetic maps that include all the aforementioned emission 
components, in the central panel of fig. || we have produced synthetic spectra extracted from 
a core (top) and an outskirts region (bottom). The core region corresponds to an angular size 
of 1° and the outskirts region is defined as a ring with inner and outer radii of 0.5° and 1.5° 
respectively (0.5° ~ 1 Mpc at the red-shift of Coma cluster). Based on the above findings, the 
idea is that by looking at the cluster core one should see the 7r°-decay and contributions 
whereas by observing the outskirts one should see IC emission from CR electrons accelerated at 




Figure 2: Left: Synthetic map of the integrated photon flux above 100 keV (top) and 100 MeV (bottom) in units 
"ph cm -2 s _1 arcmin -2 " from IC emission by secondary e^ 1 (top-left), primary e~ (top-right, bottom-right), and 
7T -decay (bottom- left). Each panel measures 15 /i _1 Mpc on a side. Center: Synthetic spectra from 10 keV up 
to 10 TeV extracted from a core (top) and an outskirts region (bottom). Right: Comparison of synthetic spectra 
with nominal sensitivity limits of future 7-ray observatories (thick-solid lines). For INTEGRAL-IBIS imagers 
(ISGRI and PICsIT) the curves correspond to a detection significance of 3a with an observing time of 10 6 s. All 
other sensitivity plots refer to a 5er significance. EGRET and GLAST sensitivities are shown for one year of all 
sky survey and Cherenkov telescopes (MAGIC, VERITAS and HESS) for 50 hour exposure on a single source. 

accretion shocks. Fig. |2| (center) shows that at high photon energy (> 100 MeV) the spectrum 
of radiation is indeed dominated by 7r°-decay (solid thin line) in the core region (top panel) 
and by IC emission (dotted line) from primary electrons in the outer region (bottom). Below 
~ 100 MeV the flux from the outskirts region is still strongly dominated by IC emission from 
primary e~ . This contribution is much reduced in the narrower field of view of the core region. 
However it is still significant and at the level of IC emission from secondary e ± . Because 
the radiation flux from these two components depends on the assumed normalizations and the 
actual shock structure subtended by the field of view, their almost equality in fig. || (center, top 
panel) is obviously only indicative. Nevertheless they are expected of comparable intensity. The 
contribution of the primary CR electrons could be estimated by measuring its radial dependence 
in the outer regions, where it dominates, and then by extrapolating its value for the core region. 

5 Discussion 

In the right panel of fig. ^ we plot sensitivity limits of planned 7-ray observatories together 
with the total radiation spectra from the core (dotted line) and outskirts (dashed line) regions. 
The sensitivity limits are plotted for point sources and will be somewhat worse for the case 
of clusters given their large angular size. In any case, according to this plot, several future 
experiments should be sensitive enough to detect the computed non-thermal emission at most 
photon energies. In particular the IBIS imager onboard INTEGRAL should readily measure the 
flux between 100 keV and several MeV. In addition GLAST and Cherenkov telescopes should 
be able to detect both the core 7-ray emission from 7r°-decay as well as the IC flux directly 
produced at accretion shocks above 100 MeV and 10 GeV respectively. 

These measurement will provide important information about GCs. First, the extended 
emission is produced at the location of of accretion shocks. Merger shocks have occasionally 
been observed in the core of clusters as relatively weak temperature jumps. However, strong 
accretion shocks have yet to be observed directly. Thus, detection and imaging of IC emission 
from primary electrons would provide an opportunity to directly observe GC accretion shocks. 

Furthermore, measuring the flux about 100 keV together with radio measurements, will 



allow an estimate of (or upper limits on) the ICM magnetic field. It is important, however, 
to separate out the contribution from the CR electrons accelerated at the outlining shocks. 
Because the magnetic field strength is expected to drop at the cluster outskirts, these electrons 
likely generate only a weak radio emission. Nonetheless, they produce a strong IC flux which 
can easily dominate the total soft 7-ray emission (cf. fig. |2|). In this case separating out their 
contribution would be an important step for correctly measuring ICM magnetic fields. 

One of the most awaited experiments is related to the measurement of the 7-ray flux at and 
above 100 MeV. This is important in order to convalidate or rule out secondary e ± models for 
radio halos in GCE3cJ'c3 and in order to estimate the non-thermal CR pressure there Q. In this 
perspective several authors estimated for nearby clusters the 7-ray flux from 7r°-decay. However, 
radiation flux from IC emission can be comparable Ej. Therefore, for a correct interpretation 
of the measurements the contributions from these two components will need to be separated as 
outlined here. Estimating the 7-ray flux from IC emission due to primary CR electrons will also 
be instrumental for addressing the contribution of this emission mechanism to the CGB&Qci 

Measurements at photon energies, from ~ 10 — 100 GeV to 1 — 10 TeV, will also be of great 
interest. Because the radiated energy spectrum is directly connected to the distribution function 
of the emitting particles, measuring the flux at different energies will provide information about 
the acceleration mechanism. In particular, the observed spectrum could differ from our predic- 
tions if the accretion shocks were modified by CR pressure. Finally the maximum energy of 
the accelerated CR electrons could be inferred from a spectral cut-off, although only as long as 
the latter turns out below the energy threshold for 77 absorption (~ TeV). Given the different 
environmental conditions with respect to supernova remnants, observations of galaxy clusters 
offer the interesting possibility for a broad exploration of shock acceleration physics. 
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